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§ ; ABSTRACT 

. We present the full source catalogue from the Australia Telescope 20 GHz (AT20G) Survey. 

The AT20G is a blind radio survey carried out at 20 GHz with the Australia Telescope Com- 
pact Array (ATCA) from 2004 to 2008, and covers the whole sky south of declination 0°. 
The AT20G source catalogue presented here is an order of magnitude larger than any previ- 
ous catalogue of high-frequency radio sources, and includes 5890 sources above a 20 GHz 
flux-density limit of 40 mJy. All AT20G sources have total intensity and polarisation mea- 
sured at 20 GHz, and most sources south of declination —15° also have near- simultaneous 
flux-density measurements at 5 and 8 GHz. A total of 1559 sources were detected in polarised 
total intensity at one or more of the three frequencies. 

The completeness of the AT20G source catalogue is 91 per cent above 100 mJy beam -1 
and 79 per cent above 50 mJybeam -1 in regions south of declination —15°. North of —15°, 
some observations of sources between 14 — 20 hr in right ascension were lost due to bad 
weather and could not be repeated, so the catalogue completeness is lower in this region. 
Each detected source was visually inspected as part of our quality control process, and so the 
reliability of the final catalogue is essentially 100 per cent. 

We detect a small but significant population of non-thermal sources that are either un- 
detected or have only weak detections in low-frequency catalogues. We introduce the term 
Ultra-Inverted Spectrum (UIS) to describe these radio sources, which have a spectral index 
a(5, 20) > +0.7 and which constitute roughly 1.2 per cent of the AT20G sample. 

The 20 GHz flux densities measured for the stro ngest AT20G source s are in excellent 
agreement with the WMAP 5 -year source catalogue of Wrig ht et al.l (2009), and we find that 
the WMAP source catalogue is close to complete for sources stronger than 1.5 Jy at 23 GHz. 

Key words: radio continuum: general - catalogues - surveys - galaxies: active - cosmic 
microwave background - methods: data analysis 



* E-mail: tara@physics.usyd.edu.au 



1 INTRODUCTION 

Large-area high-frequency radio surveys are time-consuming, and 
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Those that have been completed are either deep but covering small 
areas, or shallow all-sky surveys. As a result, our knowledge of the 
high frequency (> 10 GHz) radio source population is poor. 

In addition to the scientific benefits of studying the radio 
source population at high frequencies, large scale surveys are use- 
ful for foreground subtraction for Cosmic Microwave Background 
(CMB) experiments. Measurement of CMB anisotropics is lim- 
ited by contamination from astronomical foregrounds, both Galac- 
tic and extragalactic. To improve the efficiency of the component 
separation techniques, observations such as those of ESA's Planck 
mission are performed on a broad spectral region ranging from a 
few tenths to hundreds of GHz. At frequencies up to ~ 100 GHz, 
extragalactic radio sources are the major contaminant s on angular 
scales smaller than 30 arcminutes dde Zotti et alj|2005k so identifi- 
cation of radio sources at high frequencies is critical. 

The first blind radio survey above 8 GHz was carried out by 
iTavlor et all (1200 ll) with the Ryle telescope at 15.2 GHz. The sur- 
vey covered a 63 deg 2 region and detected 66 sources to a limiting 
flux density of 20 mJy. Waldr arrTet al.1 d2003h extended the survey 
to 520 deg 2 , detecting 465 sources to a flux density limit of 25 mJy 
(the 9C survey). Both surveys found that the existence and flux 
density of sources at 15 GHz cannot be accurately predicted by ex- 
trapolation from lower frequency radi o surveys such as t he NRAO 
VLA Sky Survey (NVSS) at 1.4 GHz dCondon et all 19981) . further 
demonstrating the need for large scale high frequency surveys for 
population characterisation and source subtraction in CMB studies. 

The Wilkinson Microwave Anisotropy Probe (WMAP) sur- 
vey, which covers t he whole sky at 23, 33, 41, 61 and 94 GHz 
(Bennett et al. 2003), is the first all-sky radio survey above 5 GHz. 
The WMAP Point Source Catalogue constru cted from the 5 yr 
maps contains 390 sources ( Wright et al. 2009), compared to 323 
and 208 so urces found in the p revious 3 yr and 1 y r maps re- 
spectively (iBennett et all 120031 : iHinshaw et"ai1 120070 . Recently, 
Mass ardiet all (2009) have detected 516 sources in the 5 yr WMAP 
maps by exploiting a combination of blind and non-blind detection 
approaches. Section[7]presents a detailed comparison of our AT20G 
results with the 5 yr WMAP point-source catalogue. 

Table Q] summarises some earlier large-area high-frequency 
radio surveys, showing the context in which the AT20G survey 
was designed. A pilot survey for the AT20G at 18.5 GHz was 
carried out in 2002 and 2003 with th e Australia Telesc ope Com- 
pact Array (ATCA) (iRicci et al.ll200i ISadler et al.ll2006h and de- 
tected 173 sou rces stronger than 100 mJy in the declination range 
-60°to -70°. lRicci et al J d2004h confirmed that ATCA (with cus- 
tom hardware) had the capability to rapidly survey the sky at high 
frequencies. The differential source counts for extragalactic sources 
from the pilot survey were found to be in good agreement with the 
Waldr am et "all d2003h 15 GHz survey. 

In this paper we present the full catalogue from the Australia 
Telescope 20 GHz Survey (AT20G). The survey covers 20,086 
square degrees (the complete Southern sky to declination 0°) to 
a limiting flux density of 40 mJybeam -1 . We followed up candi- 
date sources detected in the survey at 20 GHz and also have near- 
simultaneous follow-up observations at 5 and 8 GHz for AT20G 
sources south of declination —15°. An accompanying paper (Mas- 
sardi et al., in preparation) will provide more detailed statisti- 
cal analysis of the AT20G sample. A subset of the 320 brightest 
(S20 > 0.5 Jy) extragalac tic (| b\ > 1.5°)AT20 G sources were pre- 
sented and discussed by Massardi et al The Galactic plane 
was included in our scanning survey but no follow-up observations 
were carried out at |6| < 1. 5°, except for a blind survey of optically 
thick compact HII regions (Murphy et al. 2Q0§). 



In Section [2 we describe the survey and follow-up observa- 
tions, and in Section [3] we describe the data reduction process. In 
Section HI we calculate the accuracy of our measured positions and 
flux densities. Section \5\ presents the source catalogue and defines 
its format, with the completeness and reliability of the catalogue 
discussed in Section[6] As an additional investigation into the com- 
pleteness, Section [7] compares our catalogue with the 5 yr WMAP 
results. Finally, Section [8] discusses some statistical properties of 
the sample and Section[9]presents our conclusions. 



2 OBSERVATIONS 

The key feature of the AT20G is a two-phase observing strategy. 
The first phase of our observations exploited the fast scanning capa- 
bility of ATCA, using a wideband analogue correlator, to carry out 
a blind survey. Candidate sources from the scanning survey were 
then observed in the regular snapshot mode of the ATCA. The re- 
sults from the scanning survey, which is complete to a deeper level 
but has lower reliability, will be presented in a companion paper 
(Hancock et al., in preparation). 



2.1 Survey mode 

The first phase of our observations consisted of a raster of blind 
scans of the entire southern sky at 20 GHz, using the ATCA in 
fast scanning mode (in which it can achieve a speed of 15 degrees 
min -1 in declination at t he meridian). The A TCA has low noise 
very wideband received (Mo orev et al.ll2008h which was used to- 
gether with a custom analogue correlator. The correlator has 8 GHz 
bandwidth (Roberts, in preparation) and was originally developed 
as part o f the collabora tion for the Taiwanese CMB experiment 
AMiBA iho et alj|200lh . We used it to take dual orthogonal po- 
larization data from three of the six 22 m dishes of the ATCA. The 
lag-correlator measured 16 visibilities as a function of differential 
delay for each of the three antenna pairs used. 

Our custom correlator had no mechanism for allowing for ge- 
ometrical delay as a function of the position in the sky, so the scans 
had to be performed along the meridian corresponding to zero delay 
for the East- West configuration ATCA. All the survey observations 
used antennas 2, 3, and 4, situated on stations W102, W104, and 
W106, giving two 30.6 m (2 — 3 and 3 — 4) baselines and a single 
61.2 m (2 — 4) baseline. The shortest baseline of 30.6 m means the 
survey has reduced sensitivity to extended sources (a 50 per cent 
reduction in amplitude for source size > 45"). Larger sources will 
only be included in the survey if they have flux density in smaller 
angular size cores or hot spots above the survey limit. 

The scanning strategy consisted of sweeping sky regions 10° 
or 15° wide in declination and using Earth rotation to cover the 
full 24 hour right ascension range in a zig-zag pattern. To achieve 
Nyquist sky coverage we used multiple zig-zag scanning paths over 
multiple days - each day the scanning path was shifted by half a 
beam width. Along the scan a sample was collected every 54 ms 
(3 samples per beam), which resulted in an rms noise of 12 mJy. 
Poor quality scans (due to weather or equipment error) were re- 
peated, so that the sky coverage is as uniform as possible. The ana- 
logue correlator outputs for each set of 24 hour observations (in- 
terleaves) were combined together and calibrated to produce maps 
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Table 1. Comparison of the AT20G with other high-frequency radio surveys. Note that while the WMAP imag es cover the wh ole sky, regions at low Galactic 
latitude (|6| < 5°) are excluded from the point-source search. A reanalysis of the WMAP survey by Massa rdi et al.l ([2009) resulted in a catalogue of 516 
sources above a limit of 695 mJy. 



Survey 


Frequency 


Sky area 


Flux limit 


N sources 


Reference/s 




(GHz) 


(deg 2 ) 


(mJy) 






Ryle 


15 


60 


20 


66 


(Tavlor et al. 2001) 


9C 


15,43 


520 


25 


465 


(Waldram et al. 2003) 


9C 


15 


115 


10 




(Waldram et al. 2009) 


9C 


15 


29 


5.5 




(Waldram et al. 2009) 


WMAP 


23, 33,41,61,94 


32,177 


1000 


390 


(Bennett et al. 2003: Hinshaw et al. 2007: Wrmht et al. 2009) 


AT20G pilot 


18 


1216 


100 


126 


(Ricci et al. 2004) 


AT20G 


5, 8, 20 


20,086 


40 


5867 


(This work: Massardi et al. 2008) 



Table 2. Observation dates for the scanning observations for each of the 
declination bands. Overlapping dates are due to patching data being ob- 
served at the end of each years observations. 



| secondary A 5 min | 



Declination Band 



Observation Dates 



-30° 




> 


-40° 


11 Aug -31 Aug 2004 


-40° 


^ 5 




-50° 


20 Aug -31 Aug 2004 


-50° 


>S 




-60° 


9 Sep - 2 Oct 2005 


-60° 


^ 5 




-70° 


23 Sep - 2 Oct 2005 


-70° 






-80° 


16 Sep -20 Sep 2005 


-80° 


^ 5 




-90° 


20 & 29 Sep 2005 


-15° 






-30° 


16 Aug - 3 Sep 2006 


0° 


^ 5 




-15° 


23 Aug - 9 Sep 2007 


-85° 






-90° 


7 & 9 Sep 2007 



with an overall rms noise of ~ 10 mJy. Table[2]shows the observing 
schedule for each declination band in the scanning survey. 

The initial calibration of the interleaves was achieved by a 
daily transit observation of a nearby known calibrator. A second 
round of calibration was based on those sources from the ATCA 
calibrator catalogu^l that fell within each of the interleave obser- 
vations. Typically there were about 10 such sources within each 
24 hour observation. A third and final round of calibration was done 
using newly detected strong sources within the map itself. 

The delay steps in the analogue correlator were not precisely 
equal (as they are in a digital correlator) so the spectrum is not 
quite the same as the Fourier transform of the lags. Hence standard 
synthesis techniques using discrete Fourier transforms could not 
be used for reliable source detection. Furthermore the poor (u, v) 
coverage of the survey (30 and 60 m EW spacings only) required 
a custom source detection program that worked in a CLEAN-like 
fashion, detecting the strongest source, fitting a template dirty beam 
and recording its location and flux estimate before moving to the 
next source. Sources brighter than 5cr (~ 50 mJybeam -1 ) were 
scheduled for follow-up observations, as described in the next sec- 
tion. 



2.2 Follow-up mode 

Each of the candidate sources identified in the first phase of obser- 
vations was observed in the regular snapshot mode of the ATCA 
to confirm the detection and to measure an accurate position, flux 
density and polarisation. The flux densities from the initial scan 



block A 
20 targets 
40 sec each 



| secondary A 5 min | 
| secondary B 5 min | 



secondary A 5 min~| 



block A 
20 targets 
40 sec each 



LST +2 hr 



| secondary A 5 min | 
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Figure 1. Schematic of the follow-up observing schedule, showing blocks 
of target sources enclosed by phase calibrators. 



survey are accurate to ~ 20 per cent (see Hancock et al., in prepa- 
ration) and so we expected a fraction of the sources to fall below 
our specified survey detection threshold when re-observed. 

The follow-up observations were scheduled to be as close as 
practical (typically within a few weeks) to the initial survey obser- 
vations for a particular declination band, to reduce the effects of 
source variability on survey completeness. Table [3] lists the details 
of each follow-up run, including the shortest antenna spacing and 
resolution at each frequency. 

The 20 GHz follow-up observations were performed using a 
hybrid array configuration (i.e., one in which there were North- 
South as well as East- West baselines) with the standard ATCA dig- 
ital correlator. We used two 128 MHz bands centred at 18 752 MHz 
and 21 056 MHz and two polarisations. During data processing the 
two bands were combined to form a single 256 MHz wide band 
centered at 19 904 MHz which is the reference frequency for our 
20 GHz observations. 

Candidate sources were prioritised for follow-up observations 
in order of decreasing flux density to minimise problems with sam- 
ple completeness. The sources were scheduled in blocks of ~ 20 
targets, with a nearby secondary calibrator observed for ~ 5 min- 
utes at the beginning and end of each block. The sources within 
each block were observed in the fast mosaicing mode of the ATCA 
to reduce the slew time between pointings. We aimed to get two 
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or more 40 second cuts for each source, at different hour angles, 
so as to obtain reasonable (u, v) plane coverage. Fig. Q] shows a 
schematic of the observing schedule. Using this method up to 500 
candidate sources were observed each day. 

For declinations < — 15° each 20 GHz observing run was fol- 
lowed by lower frequency observations of the same sets of sources. 
In these runs we used an East- West extended array configuration, 
with two 128 MHz bands centred at 4 800 MHz and 8 640 MHz. 
For the rest of this paper, these frequencies will be referred to as 
'5' and '8' GHz. The lower frequency observations were conducted 
within a couple of weeks of the 20 GHz observations, so as to ob- 
tain near simultaneous spectral data for our sample. Sources in the 
most northern declination band (5 > —15°) were not observed at 5 
and 8 GHz because of the poor (u, v) coverage of EW arrays near 
the equator. 

We completed several extra observing runs, at both 20 and 
5/8 GHz for the purposes of replacing bad quality data from pre- 
vious runs. Table [3] lists the details of our observing runs, with the 
first column showing the matching high frequency-lower frequency 
epochs. The primary beam FWHM of the ATCA is 2.4, 5.5, and 
9.9 arcmin at 20, 8, and 5 GHz, respectively. 



2.3 Additional Observations 

In addition to the main survey described in this paper, we have com- 
pleted several complementary sets of observations which will be 
presented in accompanying papers. In October 2006 we observed a 
sub- sample of the bright sources to obtain high sensitivity polari- 
sation measurements (Burke- Spolaor et al. in preparation). These 
observations used the most compact configuration of the ATCA 
(H75), which provides the best coverage of a range of short spac- 
ings, and longer integration times than for the main survey, so that 
higher quality images could be made. 

Nine highly extended sources wer e selected from low fre - 
quency catalogues - PMN at 4.85 GHz dGriffith & Wrightlll993h 
and SUMSS at 843 MHz dMauch etal.l l2003) - and observed in 
mosaic mode to improve the flux density estimation at 20 GHz 
( Burke- Spolaor et "atl l2009l) . Data for seven of these sources have 
been incorporated into our catalogue in order to avoid flux density 
underestimation due to resolution effects. However, since we do 
not have equivalent 5 and 8 GHz observations of these objects we 
cannot use these data for analyses of radio spectra. 

During the AT20G follow-up survey, the 6 km antenna of the 
ATCA was operational for most of the observing time, however 
the data from the five much longer baselines could not be easily 
included in our standard processing pipeline. As a separate pro- 
gram, we estimated the 6 km visibility by taking the ratio of the 
scalar amplitude at the 6 km baseline to the scalar amplitude at 
shorter baselines. This removes the effect of atmospheric decorre- 
lation and avoids the need to phase calibrate the long baselines. 
Approximately 90 per cent of the AT20G sources have 6 km data 
available. These 6 km visibilities allows us to identify sources that 
appear point-like in the compact configuration, but extended on the 
0.3 arcsec scale corresponding to the 6 km baselines. The full anal- 
ysis and results from this high resolution follow-up will be pre- 
sented in Chhetri et al. (in preparation). An immediate application 
of this data was the identification of extended flat spectrum sources 
which are either thermal, or compact non-thermal sources extended 
by gravitational lensing. Gravitational lens candidates are being 
observed with ATCA 6 km configurations at 7 mm. The thermal 
sources identified by this procedure are Galactic objects, mainly 



planetary nebulae. They are flagged as Galactic in the main cata- 
logue - see Section [54l 



2.4 Galactic Plane Follow-up 

Neither the scanning survey nor the follow-up observations were 
well suited to imaging diffuse Galactic sources. We decided to ex- 
clude sources within |6| < 1.5° from our main follow-up survey 
and instead follow up selected subsamples of these sources with 
targeted observing runs. The first of these was the follow-up of a 
sample of ultra- and hyper-compact HII regions, selected on the 
basis of their rising spectral index between 843 MHz (from th e 
2nd Epoch Molonglo Galactic Plane Survey ( Mur phy et ai]|2007h ) 
and 20 GHz(q(0.84 3, 20) > 0.1). Further information is given in 
lMurphvetaTld2009l) . 



2.5 Variability Sample 

Many large-area radio continuum surveys have been conducted at 
frequencies of 1.4 GHz or below, where the long-term variability 
of most radio sources is low. As a result, the source catalogues 
from these surveys can continue to be used with a high level of 
confidence for many years after the survey was made. This is not 
necessarily true for surveys at frequencies above 5 GHz, where the 
source population is increasingly dominated b y flat- spectrum radi o 
sources which are expected to be variable (e.g. lCondon et al.ll989h . 

To study the variability of the 20 GHz radio- source population 
on timescales of years, we re-observed a sample of 170 sources (at 
declination —60° to —70° and with 20 GHz flux densities above 
100 mJy) at several epochs over the course of the AT20G survey. 
All the sources in this variability sample were originally d etected 
in the AT20G Pilot Survey in 2002-2003 (iRicci et al.l2004l) . giving 
up to five epochs of observation for some of these sources. 

The early ( 2002-2004) 20 GH z data for the variability sample 
are discussed by Sad ler et Id] {2006), who conclude that the general 
level of variability in sources selected at 20 GHz is relatively low on 
timescales of 1 — 2 years, with a median variability index of 6.9 per 
cent at 20 GHz over a one-year time interval. We therefore expect 
the AT20G catalogue presented here to be reasonably stable, in the 
sense that if we were to reobserve the survey area on timescales of 
a few years, would expect the new source catalogue to contain most 
of the same sources as the old one. 

The AT20G catalogue presented in this paper only lists a sin- 
gle epoch for sources in the variability sample (typically the 2004 
observation, except in cases of poor-quality data). A separate pa- 
per (Sadler et al., in preparation) will present a full catalogue and 
detailed analysis of the AT20G variability sample. 



3 DATA REDUCTION 

We developed a fully automated custom analysis pipeline to edit, 
calibrate, and image the data from the follow-up observations. This 
included a suite of quality control routines to ensure consistent data 
quality in the final catalogue. The software was developed using the 
scripting language Python, and the data reducti on was done wit h 
the aperture synthesis reduction package Miriad dSault et all 19951) . 

The data reduction was followed by a cataloguing stage in 
which the best quality sources were selected for the final cata- 
logue. This stage also included manual quality control, in which 
all sources included in the final catalogue were inspected by sev- 
eral AT20G team members (EM, EMS, JAE, RDE, TM). In the rest 
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Table 3. Follow-up observations at 20 GHz (C), to observe them at 5 and 8 GHz (O) or to repeat previous bad quality observations (R). 



Epoch 


Declination 


Freq. 1 


Freq. 2 


Array 


Shortest 


Beam l a 




Beam 2 


Dates 


Reason 




range 


MHz 


Config. 


spacing 


arcsec 




arcsec 






1 


-50° to -30° 


18752 


21056 


H214 


80 m 


10.8 x 10. 


.8 




21 Oct -27 Oct 2004 


C 


1 


-50° to -30° 


4800 


8640 


1.5C 


77 m 


8.3 x 12. 


.8 


4.6 x 7.13 


04 Nov - 08 Nov 2004 





2 


-90° to -50° 


18752 


21056 


H168 


61 m 


13.9 x 13. 


,9 




27 Oct -31 Oct 2005 


c 


2 


-90° to -50° 


4800 


8640 


1.5C 


77 m 


8.3 x 8. 


.8 


4.6 x 4.9 


12 Nov - 15 Nov 2005 





3 


-90° to -30° 


18752 


21056 


H214 


80 m 


10.8 x 10. 


,8 




29 Apr - 03 May 2006 


R 


3 


-90° to -30° 


4800 


8640 


1.5D 


107 m 


8.3 x 9. 


.5 


4.6 x 5.3 


19 Jun-23 Jun 2006 


R,0 


4 


-30° to -15° 


18752 


21056 


H214 


80 m 


2.0 x 5. 


,1 




14 Oct -17 Oct 2006 


C 


4 


-30° to -15° 


4800 


8640 


1.5B 


30 m 


8.3 x 21. 


,1 


4.6 x 11.7 


09 Nov - 12 Nov 2006 





5 


-90° to -15° 


18752 


21056 


H214 


80 m 


10.8 x 10. 


,8 




11 May -16 May 2007 


R 


5 


-90° to -15° 


4800 


8640 


1.5C 


80 m 


8.3 x 21. 


.1 


4.6 x 11.7 


04 May - 10 May 2007 


R,0 


6 


-90° to 0° 


18752 


21056 


H214 


80 m 


10.8 x 10. 


,8 




26 Oct - 30 Oct 2007 


C,R 


7 


-15° to0° 


18752 


21056 


H75 


31 m 


33.9 x 33. 


,9 




22 Aug - 26 Aug 2008 


R 



Note that the beam size changes with declination, and so these are estimates based on the typical declination of sources observed in a particular epoch. 



of this section we describe the details of the data reduction, qual- 
ity control and cataloguing. An overview of the process is given in 
Fig. El 

3.1 Flagging poor quality data 

Weather conditions can seriously affect the quality of the high- 
frequency data. Attenuation of the signal by atmospheric water 
vapour can decrease the sensitivity of the observations, and atmo- 
spheric turbulence can produce phase fluctuations that may produce 
visibility amplitude decorrelation. Hence data collected in periods 
of bad weather were removed before further processing. In partic- 
ular, calibrator data must be of high quality otherwise it introduces 
errors in the calibration solutions that affect the whole dataset. 

An atmospheric seeing monitoring system operates at the 
ATCA site. This measures the differential phase variations in a geo- 
stationary satellite signal caused by tropospheric water vapour fluc- 
tuations (as described in Middelberg et al. 2006). We used this data, 
in conjunction with the system temperature measurements from 
the antenna receivers (to estimate tropospheric opacity) to develop 
semi-automatic flagging criteria. Specifically, we discarded data 
from all the periods in which there was decorrelation greater than 
10 per cent. In cases where a calibrator was excluded, the block of 
target sources associated with that calibrator was also excluded. 

In a majority of epochs less than a few per cent of the data 
was flagged. Very occasionally, bad weather required large blocks 
of data in the follow-up survey to be edited out. Most of the time we 
were able to reobserve these blocks in clean-up runs. However, in 
the declination band —15° to 0° there are several regions between 
14 hr and 20 hr in right ascension which are still incomplete. 

3.2 Calibration 

Primary flux calibration and bandpass calibration were done in 
the standard way using PKS B 1934— 638, with the assumed fluxes 
shown in Table [4] More i nform ation about the ATCA flux scale at 
20 GHz is given in lSaultl d2003h . 

For the secondary flux calibration we followed a non-standard 
procedure which is summarised in Fig. [3] In each epoch of our 
observations we typically observed around ~ 50 secondary cali- 
brators. To calculate an accurate flux density for each secondary 
calibrator we calculated the mean of the individual snapshot flux 



Table 4. Assumed fluxes for the primary calibrator PKS B 1934-638. 



Frequency 


Flux 


(MHz) 


(Jy) 


4800 


5.83 


8640 


2.84 


18752 


1.04 


21056 


0.88 



densities across the whole run, excluding snapshots which had a 
flux density greater than two standard deviations away from the 
mean. Each target source was then calibrated using the secondary 
calibrator associated with its observing block. 

Our observation strategy meant that we had insufficient data to 
determine the instrumental polarisation corrections from the sec- 
ondary calibrators. Hence we calculated the leakage terms using 
the primary calibrator, PKS B 1934—638. The linear polarisation of 
this calibrator is known to be not variable and less than 0.2 per cent 
of the total source flux density at each of our observing frequen- 
cies. To determine the leakage terms it was assumed to be unpo- 
larised. We adopted these leakage values for all the secondary cal- 
ibrators, simultaneously calculating the time-dependent complex 
antenna gains, the residual xy phase differences (x and y are the 
orthogonal linear polarisations), and the Q and U Stokes parame- 
ters of the calibrators. The polarisation calibration was then applied 
to the target sources and /, Q, U and V Stokes parameters were de- 
termined for all of the target sources. 

3.3 Flux density measurements 

3. 3. 1 Triple product fluxes 

We measured flux densities at all three frequencies using the triple 
product method implemented in the Miriad task CALRED. The am- 
plitude of triple product is the geometric average of the visibility 
amplitudes in a baseline closure triangle 

A T p = \A4i, 2 • A 2 ,3 -A 3il (1) 

and its phase is the closure phase. Compared to measuring the flux 
densities from images, as is typical, this method of measuring flux 
densities is robust to the effects of phase decorrelation (Ricci et 
al., in preparation). This meant that we were able to recover flux 
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Figure 2. An overview of the AT20G custom data reduction process. See 
Fig[3]for more detail about calibration. 



density measurements for 143 of our target sources that could not 
be imaged due to poor weather conditions (these are marked as 
'poor' in the catalogue). 
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Figure 3. The custom calibration process developed for the follow-up data 
reduction. 



3. 3. 2 Extended source fluxes 

Sources that are partially resolved, or have multiple components, 
are not well-characterised by the triple product flux measurement. 
We identified extended sources in our sample using several criteria 
based on the ratio of the triple product flux density to the flux den- 
sity measured on the shortest baseline. We visually inspected the 
20 GHz visibilities and images for each of these sources, to con- 
firm a sample of 337 extended sources. Fig.|4]shows an example of 
a typical AT20G extended source. 

For the extended sources, we calculated the integrated flux 
densities at 5, 8 and 20 GHz using the amplitude of the signal mea- 
sured by the shortest baseline. Any source extended at 20 GHz was 
assumed to be extended at 5 and 8 GHz. In most cases the shortest 
physical baseline used in the follow-up observations is either 60 m 
or 80 m (see Table [3]), so our 20 GHz flux densities for sources 
larger than 20 arcsec will still be underestimated using this method. 
Although improved algorithms could be used, the snapshot obser- 
vations are inadequate to make reliable estimates. The extended 
sources will be discussed in more detail in a future paper. Note that 
the H75 array observations used for epoch 7 have closest spacing 
of 30 m so the short spacing flux estimate is good for sources up to 
about 1 arcmin. 

For extended sources with multiple components larger than 
40 arcsec, the shortest spacing can correspond to a minimum in the 
visibility and is not a useful estimate of total flux. For these sources 
we quote the triple product flux density that would correspond to 



the flux density of the dominant component and flag the flux as 
'poor' in the source catalogue. 

Nine of the extremely extended sources were observed 
separately and we have u sed the flux densities determined in 
iBurke-Spolaor et aD d2009h for seven of them. These are discussed 
in Section [57TI 

3.4 Imaging 

We imaged all sources using the standard Miriad process, and de- 
convolved them with a small number of CLEAN iterations (typi- 
cally 50). This was done primarily for the purpose of visual inspec- 
tion, since all flux densities were determined using the triple prod- 
uct as discussed in the previous section. Fig. |4] shows examples of 
a typical point source, AT20G J00444 1-353034, and a typical ex- 
tended source, AT20G Jl 10622— 210858, at all three frequencies 
plus the total polarized intensity. 

3.5 Polarisation 

Polarization maps in the four Stokes parameters were created using 
the Miriad routine IMPOL. The total polarized intensity (P) for each 
source was measured by calculating the fractional polarization (m) 
in the image at the position of the source, and multiplying it by the 
flux density calculated using the triple product. The error on the 
polarized flux was calculated using 

P err = V2— (2) 
X 

where a n is the error associated with the noise, which was calcu- 
lated from the rms noise measured in the Stokes V image, and % is 
ratio of the flux density in the Stokes I image to the flux density as 
calculated from the triple product. This corrects for phase decorre- 
lation, which would affect both polarized and total flux density in 
the same way. 

The polarised flux was measured at the position of the peak 
in the Stokes I image. This will be equivalent to the integrated po- 
larization for the unresolved sources but will be neither the peak 
nor the integrated polarization for the extended sources. For the 
seven extremely extended sources discussed in Section [5Jl w e have 
used t he integrated polarization determined by Bur ke- Spolaor et all 
b009h . 

We used the sources observed multiple times as part of the 
variability sample to investigate the reproducibility of our polarisa- 
tion measurements (see Section [87Tb . Based on these experiments 
we developed the following rules for defining a detection 

• if P ^ 3P e rr then P = 3P err as a limit 

• if m < 1% then P = O.OIS20 as a limit 

• if P < 6 mJy then P = 6 mJy as a limit 

where P is the total polarized intensity and m is the fractional po- 
larization. If no polarization detection was made, then the limit on 
the polarized flux density is the maximum of these three limits. 

3.6 Cataloguing 

All good quality observations were selected for the final source cat- 
alogue. In cases where a source was observed in multiple epochs 
(usually due to there being a poor quality observation in the ear- 
lier epoch) we selected the best source based on the data quality, 
the presence of near- simultaneous 5/8 GHz data and the offset of 
the source from the imaging pointing centre. For sources that were 
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Figure 5. Image of AT20G J00444 1-353034 (see Fig.gJ, created from the 
scanning survey data. 

observed multiple times as part of the variability study, we selected 
the first epoch of observation except in cases where that observation 
was of poor quality. 

Due to the poor (u, v) coverage of the scanning survey, a 
small percentage of the sources scheduled for follow-up observa- 
tions were in fact sidelobes of the main source. Fig. [5] shows a 
typical AT20G source as detected in the images created from the 
scanning survey. The scanning survey has only two different base- 
lines (EW 30 m and 60 m) so the images have very high sidelobes. 
In some cases bright sidelobes were inadvertently followed-up. 
We excluded sidelobes automatically by grouping all 'duplicate' 
sources within a radius of 100 arcsec (the effective beam size of the 
scanning survey) and filtering to keep only the strongest source in 
the group. Cross-matching with low frequency surveys also helped 
to confirm these selections. 

All sources were then inspected by several members of the 
AT20G team, using an online annotation tool. We used this as a 
final quality control procedure, and to identify any remaining side- 
lobes. We also used this process to identify extended sources that 
were not classified by our automatic selection criteria. 



4 ACCURACY 

4.1 Position uncertainties 

The errors in position (cr a , as) were calculated using 

°a — a a,ca\ + a a,n (3) 
&6 = Cr|cal + 0"l,n ( 4 ) 

where cr ca i is the error associated with calibration, and <j n is the 
error associated with the noise. 

The calibration term was estimated by comparison with the the 
Inter national Celesti al Reference Frame (ICRF) defining calibra- 
tors dFev et al.| [2004). The VLBI-measured positions in the ICRF 
catalogue are accurate to a milliarcsecond or below, so any dis- 
crepancy between the positions of our target sources and ICRF 
positions can be assumed to be due to calibration positional er- 
rors in our sample. There are 25 1 ICRF calibrators in our sample, 
and the positional offsets for these are shown in Fig. [6] The mean 
offsets in right ascension and declination are (Act) = 0.1 arcsec 
and (AS) = 0.0 arcsec, showing we are in good agreement with 
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Figure 4. Follow-up images for a typical AT20G point source, J00444 1—353034 (top) and a typical extended source, Jl 10622— 210858 (bottom). From left 
to right they show the total intensity at 20 GHz, 8 GHz and 5 GHz, as well as the total polarization map. Note that each image has a different intensity scale. 
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Figure 6. Offset of AT20G 20 GHz positions from ICRF positions for the 
25 1 ICRF calibrators in our catalogue. 



the ICRF positions. The values for rms scatter which we used as 
the calibration terms in Equation [3] are cr a , C ai = 0.8 arcsec and 
o~5,ca\ = 0.9 arcsec. 

The calibration term dominates over the noise for our data, are 
we assumed a common calibration term for all the target sources. 
From this we calculated the mean positional errors in right ascen- 
sion and declination for the full sample, giving cr a = 0.9 arcsec 
and as = 1.0 arcsec. Note that this method of calculating posi- 
tional errors may underestimate the error for extended sources. 



4.2 Flux density uncertainties 

We calculated the errors in the flux density measurements by 
adding in quadrature the error associated with calibration (gain er- 



The gain error is a multiplicative term (i.e., it is proportional 
to the source flux density) and is a measure of the gain stability 
over time. We estimated cr ga in for each observational epoch and 
frequency from the scatter in the visibility amplitudes of the cali- 
brators in each observing run. The mean values for the gain errors 
were found to be of order a few per cent. The noise term is an addi- 
tive term related to the interferometer noise, which is proportional 
to the system temperature. Since no source has significant flux in 
Stokes V, the rms noise levels in the V images have no gain error 
and were used as an estimate of a n . The error in the flux density 
for each source is given in the main catalogue (Table 4). Typically 
the error is 4 — 5 per cent of the total flux density. 

For extended sources the error was multiplied by the square 
root of the number of baselines ribase (normally 10 for our 5- 
antenna follow-up arrays) to correct for the fact that the flux den- 
sities for these sources are estimated using only one (the shortest) 
baseline instead of ribase- 

Many of the sources in the Australia Telescope calibrator cat- 
alogue were observed as target sources as part of the AT20G pro- 
gram. We used these sources to check the integrity of our observing 
and data reduction procedures. There are 362 AT calibrators in the 
final AT20G catalogue. For each of these we extracted the AT cali- 
brator catalogue flux density measurement that was closest in time 
to our AT20G observations of that source. Fig.|7]shows the AT cali- 
brator catalogue to AT20G flux density ratio for these sources. The 
median AT20G/AT flux density ratio is 1.06 ± 0.03. All of the out- 
liers were found to be either extended or highly variable sources. 

A small number of our target sources are included in the VLA 
Calibrator Manual We used these as an additional check on our 
flux density measurements, particularly for sources near the equa- 
tor. The VLA calibrator data is only available up to the year 2000. 
This means there is a significance difference in the times that these 
measurements were made. The scatter we see in our sample is 



' |http : / /www . via . nrao . edu/ astro/ calib /ma nual] 
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Figure 7. The ratio of AT20G 20 GHz flux density to Australia Telescope 
Calibrator Catalogue flux density, for the 362 AT calibrators in our cata- 
logue. The AT calibrator flux density closest in time to our AT20G obser- 
vations was chosen. 



quite high, due to variability of the sources over this 8 year pe- 
riod. The median AT20G/VLA flux density ratio for these sources 
is 1.08±0.10. 

We have explored the reasons for the possible flux scale differ- 
ences between AT20G, and ATCA and VLA calibrators. However, 
without simultaneous measurements on both telescopes it is diffi- 
cult to obtained a conclusive result. We see no evidence that our 
catalogue flux densities are not on the scale set by the assumed 
flux of 1934—638. We note that our fluxes are consistent with the 
WMAP flux scale (see Section[7} and strongly recommend a future 
program of simultaneous measurements of calibrator sources to tie 
the flux scales in the north and south hemispheres, and between the 
major radio observatories. 



(17-18) Polarized intensity at 8 GHz and error in mJy. 

(19) Fractional polarisation at 8 GHz. 

(20) Position angle at 8 GHz. 

(21-22) Polarized intensity at 5 GHz and error in mJy. 

(23) Fractional polarisation at 5 GHz. 

(24) Position angle at 5 GHz. 

Column 11 contains a quality flag, either g (good, 5501 sources) 
or p (poor, 389 sources). A poor quality flag indicates that there 
was lower quality data in that observation, or that, in the case of 
an extended source, the triple product flux was used rather than the 
shortest baseline flux. Hence the flux density measurement for poor 
quality sources may not be as reliable. 

Column 12 contains other flags identifying the following 
source types or issues: 

e Source is extended and the shortest baseline flux has been used 
(see Section [Oi l. 

h Source identified as a Galactic HII region (Section \5A[ . 
p Source identified as a Galactic Planetary Nebula (Section [5~4t . 
m Source identified as part of the Magellanic Clouds (Sec- 
tionlPl. 

1 Source has no match in the low frequency surveys (NVSS and 
SUMSS) (SectionES. 

b Source is large and extended (see Table |5}. The data 
in the catalogue comes from the observations presented in 
Burk e- Spolaor et all J2009I) and discussed in Section [57TI 

AT20G catalog ue sources sho uld be referred to by their full 
IAU designation (lLortet et al .11 19941) . These are of the form AT20G 
JHHMMSS-DDMMSS where AT20G is the survey acronym, J 
specifies J2000.0 coordinate equinox, HHMMSS are the hours, 
minutes and truncated seconds of right ascension, — is the sign of 
declination and DDMMSS are the degrees, minutes and truncated 
seconds of declination. 



5 SOURCE CATALOGUE 

The primary AT20G source catalogue gives the flux density and 
polarisation measurements at 20, 8 and 5 GHz, as well as the epoch 
of observation and quality flags. Table 7 shows the first 30 sources 
in the catalogue, the full catalogue is included in the Supporting 
Information for this paper, and available online through VizieQ 
The columns are: 

(I) AT20G source name. 

(2-3) Right ascension and declination (J2000). The mean errors 
(a a = 0.9 arcsec, as = 1.0 arcsec) are derived in Section|4] 

(4-5) Flux density at 20 GHz and error in mJy. 

(6-7) Flux density at 8 GHz and error in mJy. 

(8-9) Flux density at 5 GHz and error in mJy. 

(10) Epoch of observation for the three frequencies (20, 8 and 
5 GHz). 

The epoch codes are listed in column 1 of Table [3] 

(II) Quality flag (see description below) 
(12) Other flags (see description below) 

(13-14) Polarized intensity at 20 GHz and error in mJy. 

(15) Fractional polarisation at 20 GHz. 

(16) Position angle at 20 GHz. 

4 http://vizier.u-strasbg.fr 



5.1 Bright extended sources 

As discussed in iMassardi et al.l d2008h . a small number of highly 
extended sources were expected to have a 20 GHz flux density 
above our bright source sample cutoff of 0.5 Jy, but were unde- 
tected or had diminished levels of observed emission in the AT20G 
follow-up observations due to the sources' extent beyond our ob- 
serving resolution and field of view (~ 2.4 arcminutes). These 
sources were identified using extrap olated flux densities from the 
5 GHz PMN ([Griffith & Wright! 19931) and the 0.843 MHz SUMSS 
jMauch et alJl2003h data. Integrated flux density and polarisation 
measurements for these sources (shown in Table [5} were measured 
from mosaic observations done during the 2006 October polarisa- 
tion follow-up run. Their flux densities have been included in our 
main catalogue and are flagge d with a b. The observation s and data 
analysis are discussed fully in iBurke- Spolaor et al.l (120091) . The ob- 
jects with mosaiced measurements are limited to sources south of 
5 = -30°. 

We note that as discussed in Section 12.11 our survey will be 
incomplete for sources larger than 45 "and based on simple source 
count arguments many hundred sources larger than a few arc min- 
utes will be missing. Since non-thermal extended sources have 
steep radio spectra we suggest that low frequency catalogues such 
as NVSS or SUMSS rather than AT20G be used if complete sam- 
ples of the extended sources are important. 
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Table 5. Properties of 7 bright extended sources from Burke-Spolaor et al.l 
( 2009) that have been included in the AT20G source catalogue. Core and 
total integrated flux density (or limits) and errors are given in Jy. 



AT20G Name 


Source Name 


Score (Jy) 


Sis (Jy) 


AT20G J013357- 


-362935 


PKS 0131-36 


0.03±0.01 


>0.44 


AT20G J051949- 


-454643 


Pictor A 


1.32±0.04 


6.32±0.11 


AT20GJ132527- 


-430104 


Centaurus A 


5.98±0.17 


>28.35 


AT20GJ134649- 


-602430 


Centaurus B 


5.02±0.06 


8.89±0.43 


AT20GJ161505- 


-605427 


PKS 1610-60 


0.14±0.05 


2.11±0.04 


AT20G J215706- 


-694123 


PKS 2153-69 




3.40±0.21 


AT20G J235904- 


-605503 


PKS 2356-61 


0.09±0.03 


1.64±0.05 



5.2 Comparison with bright source sample 

Some bright sources (S20 > 0.5 Jy) in this catalogue have slightly 
different parame ters than were given in the Bright Source Sam- 
ple (BSS) paper (Massar di et al.l2008h . There are 23 sources which 
have a flux density that differs by more than 10 per cent from the 
BSS results. In most cases this is due to a different observational 
epoch being selected for the source in the BSS than for in the cur- 
rent catalogue. The largest group of sources in this category are 
sources from the variability sample discussed in Section |231 In the 
selection of sources for this catalogue we gave preference where 
possible to the epoch with simultaneous 5 and 8 GHz observations. 
However, in the BSS the best quality epoch at 20 GHz was chosen 
regardless of the availability of follow-up at other frequencies. 

5.3 Magellanic Clouds 

The AT20G survey area includes the Large and Small Magellanic 
Clouds, and we would expect to find some radio sources in these 
galaxies, which are effectively foreground sources for cosmolog- 
ical studies. In e arlie r work using the A T20G pilot survey data, 
iRicci etal.lJ2004h and lSadler et~aiT(l2006h simply removed 5° x 5° 
regions of sky around the LMC and SMC to avoid the problem 
of foreground contamination. However, since the compact, high- 
frequency LMC/SMC radio sources are likely to be interesting in 
their own right, we have now attempted to identify them individu- 
ally in the final AT20G catalogue. 

We searched two areas of sky defined by the following 
(J2000) coordinates: 

SMC 00:30 < a < 01:30 -71°< 8 < -75° 
LMC 04:45 < a < 06:00 -66° < 8 < -72° 

There are six AT20G sources in the SMC region and 21 in the 
LMC region. As with Galactic sources, some LMC and SMC radio 
sources may have complex, extended emission which is not well 
imaged by the AT20G snapshot observations. In such cases, the flux 
densities and positions listed in the AT20G catalogue should be re- 
garded as no more than indicative. Hence one of the main reasons 
for identifying LMC/SMC sources is to have the option of exclud- 
ing them from later analysis of the extragalactic AT20G sample. 

To determine whether each 20 GHz source was likely to be 
associated with the LMC or SMC we searched the NASA Ex- 
tragalactic Database (NEDjfl for known LMC/SMC objects at or 
near the AT20G position, and cross-matched the AT20G posi- 
tions with lower-frequency radio catalogues from iFilipovic et all 

5 |http : / / nedwww . ipac . caltech . edu| 



Table 6. AT20G sources which are probably associated with objects in the 
Large and Small Magellanic Clouds. HII = HII region, PWN = pulsar wind 
nebula. Objects in the LMC and SMC regions which we classified as back- 
ground sources are also listed. 



AT20G name AT20G position (J2000) Notes 



SMC objects 



J012407- 


-730904 


01 9zL 07 Q9 


7^ OQ OzL 1 
— id yjy yj^t. 1 


TR A ^01998 7^9zL HTT 


J012930- 


-733311 


01 9Q ^O 07 


l^K X\ 1 1 ^ 

— ID DD ll.J 


TP A <5 01 98^ 7^ZLQ T-TTT 
iiv/\o uizoj — i D^y, nn 


LMC objects 








J045153- 


-692329 


OzL S1 S^ 9S 


f.Q oa 90 A 
— \jy LD Ly 


TPA<5 0zlS91 6Q98 T-TTT 
llv/\0 LH-JZl — O7Z0, nil 


J050950- 


-685305 


OS OQ SO f\\ 
yjj yjy ju.ui 


68 S^ OS f\ 
— Do JD UJ.O 


Part nf TsJHP 1 8S8 T-TTT 
r al I Ul rNvJv.^ 10J0, nil 


J051317- 


-692222 


OS 1 ^ 1 7 7S 
U J ID 1 / . / J 


AO 99 99 f. 
— \jy zz zz.o 


Large HII region complex 


J052212- 


-675832 


OS 99 19 f\l 
yjj zz iz.u / 


fH SR X) Q 
— O /Jo DA.y 


Part nf TsJOP 1 T-TTT 
r dlL Ul 1MUU iyD\J, nn 


J053745- 


-691010 


OS XI zLS S1 

UJ D 1 HO . D 1 


— f\Q 10 10 
— yjy 1U lu.U 


P^R T0S^7— f\Q PWN 

rjlv jyjDD I — U", r VV IN 


J053845- 


-690503 


OS XA zLS f\f\ 

UJ DO tJ.DD 


f\Q OS O^ 1 
— yjy yjj \JD. 1 


Part nf TsjnP 9070 T-TTT 
rd.lL Ul 1NUU zu/u, nn 


J053937- 


-694526 


OS XI 41 

UJ Dy D 1 . L r 1 


— f\Q 4S lf\ A 
— yjy L rD zu.h- 


Part nf NOr 907Q T-TTT 

rail Ul INvJv^ ZU/", nil 


J053945- 


-693839 


uj jy 4 j. j / 


— yjy Do Dy.Z 


JNLrL. zUoU, nil 


J054004- 


-694438 


05 40 04.76 


-69 44 38.6 


Part of NGC 2079, HII 


J054011- 


-691953 


05 40 11.09 


-69 19 53.4 


PSRB0540— 69.3, PWN 


J054024- 


-694014 


05 40 24.69 


-69 40 14.5 


IC 2145, HII 


Background sources 








J004047- 


-714559 


00 40 47.90 


-71 45 59.6 




J005611- 


-710707 


00 56 11.34 


-71 07 07.0 




J011049- 


-731428 


01 10 49.61 


-73 14 28.2 




J011132- 


-730209 


01 11 32.25 


-73 02 09.9 




J045551- 


-690209 


04 55 51.56 


-69 02 09.5 




J045608- 


-701433 


04 56 08.67 


-70 14 33.3 




J050551- 


-695116 


05 05 51.89 


-69 51 16.5 




J051129- 


-680618 


05 11 29.40 


-68 06 18.1 




J051222- 


-673220 


05 12 22.54 


-67 32 20.5 




J051537- 


-672128 


05 15 37.36 


-67 21 28.4 




J051832- 


-693520 


05 18 32.48 


-69 35 20.6 




J052635- 


-674909 


05 26 35.05 


-67 49 09.2 




J054317- 


-662655 


05 43 17.54 


-66 26 55.8 




J054750- 


-672801 


05 47 50.02 


-67 28 01.8 





( 2005b and lPavne etai] (E004). who independently classified many 
of these sources as either LMC/SMC or background objects. We 
also checked for optical nebulosity of the sources by inspecting 
overlays of AT20G contours on optical images from SuperCOS- 
MOS. 

In most cases, it was possible to distinguish foreground 
LMC/SMC objects from background AGN with a high level of 
confidence. We classified 14 of the 27 sources as background extra- 
galactic objects and 13 as LMC/SMC radio sources (listed in Table 
[6]and identified in the catalogue with a flag 'ra'). As expected most 
of the LMC/SMC objects are HII regions, but the AT20G catalogue 
also includes two LMC pulsar wind nebulae, PSR J0537— 69 and 
PSRB0540-69.3. 



5.4 Galactic sources 

The Galactic plane \b\ < 1.5° was excluded from the main AT20G 
follow-up. However, as part of our analysis some sources at higher 
Galactic latitude have been identified as Galactic Planetary Nebulae 
or HII regions. These have been identified in the source catalogue 
(with flags p and h respectively) to allow them to be excluded from 
extragalactic studies. 

The 65 sources identified as PNe were found by looking for 
flat spectrum sources that were slightly extended, in terms of the 
6 km visibility data. A full analysis including this data will be 
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published separately. The 6 HII regions were identified through a 
search of the literature. 



6 COMPLETENESS AND RELIABILITY 

The overall completeness of the AT20G catalogue is a function of 
both the completeness of the original scanning survey, and the com- 
pleteness of the follow-up survey. The methods for estimating com- 
pleteness for both surveys are discuss ed in this sect i on. In Section|7] 
we compare our catalogue with the Wright et al ] d2QQ9h WMAP 
catalogue to assess completeness and reliability. 

6.1 Scanning survey completeness 

In order to estimate the completeness of the survey catalogue (the 
catalogue of candidate sources that was used as a basis for the 
follow-up observations) false point sources with known bright- 
ness and positions were inserted into the time ordered uncalibrated 
scans. The false sources were injected by taking the primary cal- 
ibrator observation, scaling it to the correct flux and declination, 
applying an inverse calibration appropriate to the time of the obser- 
vation, and adding this to the raw data. The data was then processed 
using the same machinery that created the initial sky maps and a 
list of candidate sources including the injected sources were then 
detected. Finally, the extracted sources were compared to the list 
of input sources to measure the fraction of sources recovered as a 
function of flux density. Further details of this process are discussed 
in Hancock et al. (in preparation). 

This process allowed us to assess the completeness for point 
sources. The AT20G survey is not sensitive to sources with angu- 
lar sizes much larger than 45" and so our catalogue will be in- 
complete for extremely extended sources. This brightness sensi- 
tivity limit means the AT20G survey can not detect the emission 
from most nearby spiral galaxies, even those with integrated flux 
densities above lOOmJy at 20 GHz. The nearby spiral NGC253 
is detected as an extended AT20G source (J004733-251717) with 
a 20 GHz flux density of 608 mJy, and in this case the observed 
emission appe ars to arise from a central starbur st rather than 
an AGN fsee iTingavl I2OO4I : iBrunthaler et al.ll2009h . At z=0.0008 
(d=3.3 Mpc), NGC 253 is the lowest-redshift galaxy detected in the 
AT20G survey. Although several other spiral galaxies are detected, 
including NGC 1068 (J024240- 000046; z=0.0038), NGC 4594 
(J123959-113721; z=0.0036) and NGC 4945 (J130527-492804; 
z=0.0019), the detected radio emission in these galaxies appears to 
be associated mainly with an active nucleus rather than processes 
related to star formation. 

6.2 Follow-up survey completeness 

The completeness of the follow-up survey is a function of the com- 
pleteness of the scanning survey catalogue, and of the number of 
objects that were scheduled and observed in follow-up mode. The 
number of sources that we were able to observe was limited by 
the finite length of the project and the weather conditions during 
each observing run. The brightest sources in each region of the sky 
were observed first so as to maximize the number of sources con- 
firmed and thus conserve the completeness of the follow-up obser- 
vations. Bad weather and occasional hardware malfunctions meant 
that even though a source was scheduled for observation no good 
data was obtained. Much of this missing data was able to be re- 
couped at a later stage via clean up observations, but not all. 



Fig. [8] shows the source distribution for the AT20G follow- 
up survey source catalogue. We calculated the follow-up survey 
completeness (Cf) using 

C f = ^ (6) 

where C s is the completeness of the scanning survey (estimated us- 
ing the source injection method described above); nj is the number 
of sources detected in the follow-up survey; and n s is the number of 
real sources expected from the scanning survey. Note that the flux 
densities determined from the original scanning survey were accu- 
rate to rsj 20 per cent. Hence towards the flux density cutoff there 
is a reduction in completeness due to some sources being omitted 
from the original candidate lists. 

For regions south of 8 = —15° we estimate the survey com- 
pleteness to be 91 per cent above 100 mJy and 79 per cent above 
50 mJy. The region — 15° < S < 0° is notably less complete than 
the rest of the survey. This region was observed last, and was ham- 
pered by bad weather. We had several catch-up runs but this was not 
enough to fill in the missing area within the time limits we had set 
for near- simultaneous follow-up. The main region this affected was 
between 14 hr and 20 hr in right ascension. A full analysis of the 
completeness statistics, including values for each declination band 
with ^ 2 hr right ascension zones will be presented in Massardi et 
al. (in preparation). 

Since every object detected in the follow-up survey passed ex- 
tensive quality control criteria, in addition to visual inspection, the 
reliability of the catalogue is essentially 100 per cent. 



7 COMPARISON WITH THE WMAP 5-YEAR SOURCE 
CATALOGUE 

Foreground point- source removal is an important step in the analy- 
sis of CMB anisotropy data, since the strongest sources need to be 
identified and masked out of the CMB images. The WMAP team 
have therefore produced a catalogue of bright point sources in the 
WMAP sky maps at 23, 33, 41, 61 and 94 GHz ilWright et al.l2009h 
so that these can be masked out for CMB analysis. 

Since the AT20G survey overlaps in frequency with the lowest 
WMAP band (20-25 GHz), and the AT20G has significantly bet- 
ter sensitivity and resolution than the WMAP images, the AT20G 
catalogue provides an independent check of the completeness and 
reliabi lity of the WMAP p oint- source catalogues. 

Massardi et al. ( 2009) have recently published a detailed anal- 
ysis which uses data from the AT20G Bright Source Sample 
( Massardi et al to evaluate several different strategies for 

foreground source detection in the WMAP 5 -year maps. The re- 
lease of the full AT20G catalogue will allow similar studies to push 
below the 0.5 Jy BSS limit and provide deeper complete samples 
for Planck and other CMB experiments. 

A detailed comparison of the AT20G and WMAP source cat- 
alogues is beyond the scope of this paper, so we focus here on 
two questions: (i) what is the completeness and reliability of the 
WMAP source catalogue? and (ii) how consistent are the AT20G 
and WMAP flux density scales at 20-25 GHz? 

7.1 Matching the AT20G and WMAP source catalogues 

The WMAP 5 -year source catalogue (Wright etal. 2009) con- 
tains 390 sources, of which 186 are in the southern hemisphere. 
Wrig ht et all d2009h note that the WMAP catalogue is expected to 
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be complete for sources stronger than 2 Jy in regions of the sky 
away from the Galactic plane, but also contains some sources with 
flux densities as low as 0.5 Jy at 23 GHz (WMAP K-bandfl 

In matching the final AT20G catalogue with the WMAP 5- 
year source c atalogue, we adopted the same 21.35 arcmin cutoff 
radius used bv lMassardi et all d2008h for the AT20G Bright Source 
Sample. 180 of the 186 southern sources in the WMAP cata- 
logue were matched with at least one AT20G source within this 
21.35 arcmin radius. Based on our determination of the surface den- 
sity of bright AT20G sources, we expect all these matches to be 
genuine associations. 

In 16 cases (i.e. 9 per cent of the WMAP source catalogue), 
two or more AT20G sources make a significant (> 10 per cent) 
contribution to the total flux density in the 0.93° WMAP beam at 
23 GHz. Some of these correspond to AT20G detections of several 
components (e.g. core and hotspot, or two hotspots) of a single ex- 
tended radio galaxy, while others appear to be unrelated pairs of 
sources. 



6 IWright etakl d2Q09l) note that sources with flux densities below 1 Jy are 
unlikely to be detected by WMAP unless they have 'benefited' from a pos- 
itive noise or CMB fluctuation. This in turn leads to a bias in the WMAP 
catalogue at low flux densities (see, e.g. Eddington 1913; Jauncev 1968). 



Fig. |9] shows the offsets between the AT20G and WMAP po- 
sitions for sources in common. In general these are consistent with 
the quoted position errors (typically ^4 arcmin for WMAP and less 
than 1 arcsec for AT20G). As expected, the offsets are generally 
smaller for stronger WMAP sources. 

7.2 Completeness and reliability of the WMAP 5-year source 
catalogue 

A simple comparison of surface densities shows that the WMAP 
and AT20G catalogues have similar levels of completeness for 
strong (>1.0Jy) sources at 20-25 GHz. The surface density of 
sources stronger than 1.0 Jy is 28.5±2.2sr _1 for AT20G and 
26.8=bl.7sr _1 for the WMAP K-band, where the quoted errors 
are set by the sample size in each cas^ll Since the WMAP sur- 



7 The 5-year WMAP catalogue covers 78 per cent of the sky (the remain- 
ing 22 per cent, at low Galactic latitude, is masked out) and lists 262 sources 
with 23 GHz flux densities ^1.0 Jy beam -1 and 123 sources with flux den- 
sity ^1.5 Jy beam -1 . The AT20G catalogue covers 48 per cent of the sky 
(declination < 0° with the region \b\ < 1.5° masked out) and lists 174 
sources with 20 GHz flux densities ^1.0 Jy beam -1 and 85 sources with 
flux density ^ 1 .5 Jy beam - 1 . 



AT20G Source Catalogue 13 



(DDO GBD 

oo ooo o 

CS® O GMD OTXDO O 

oo co an> 
-o © oo o oo 
-co oo • o o o 

-CD O O O O 

O O GO 



5 10 15 20 

AT2QG-WMAP position offset (arcmin) 



^ 10 



a 
o 

cv 

H 
< 



X 



C? O / (§, 




o. oo 



1 10 
WMAP5 23 GHz flux density (Jy) 



Figure 9. Offset between AT20G and WMAP positions for southern sources 
in the 5 -year WMAP source catalogue ( Wrigh t et alJl2009h . Open circles 
show WMAP sources that were matched to a single AT20G source, and 
filled circles objects in which more than one AT20G source makes a signif- 
icant (> 10 per cent) contribution to the flux density in the WMAP beam 
at 23 GHz. 



face density at 23 GHz is 94 per cent of the AT20G value, and 
the AT20G catalogue is essentially complete for point sources 
strong than 1 Jy, we estimate that the WMAP catalogue is roughly 
94 per cent complete above 1 Jy at 23 GHz. Excluding the area 
masked in the WMAP source catalogue, we find only one strong 
(> 1.5 Jy) AT20G source which is not listed in the WMAP cata- 
logue, AT20G J142432-491349. 

Appendix A lists the eight sources in the WMAP catalogue 
which have no AT20G match stronger than 250 mJy within 21.35 
arcmin of the WMAP position. Three of these (the radio galaxy 
Fornax A, a Galactic HII region and the Galactic planetary neb- 
ula NGC 7293) are highly-extended sources which are resolved out 
by the 2 arcmin ATCA beam, and so represent an incompleteness 
in the AT20G sample. Two sources are found in both surveys but 
with a very large position offset, and three sources catalogued by 
WMAP are not confirmed by AT20G and appear to be spurious. In 
general, however, there is very good agreement between the AT20G 
and WMAP source catalogues. Since the AT20G catalogue is es- 
sentially 100% reliable, and -98% of 23 GHz WMAP sources are 
matched in AT20G, we conclude that the 23 GHz WMAP catalogue 
is -98% reliable. 



7.3 Comparison of the AT20G and WMAP flux density scales 

Fig. [lOl compares the catalogued AT20G (20 GHz) and WMAP 
(23 GHz) flux densities for the 1 80 sources in common. 

To test for consistency of the AT20G and WMAP flux-density 
scales, we compared the flux densities of the 119 sources which 
were stronger than 1.0 Jy in the WMAP catalogue (to minimize 
measurement errors) and had only a single AT20G counterpart (to 
exclude extended sources for which the AT20G flux densities may 



Figure 10. Comparison of AT20G (20 GHz) and WMAP (23 GHz) flux 
densities. As in Fig. [9] open circles represent WMAP sources matched to 
a single AT20G source and filled circles objects with two or more AT20G 
sources in the WMAP beam. The horizontal dashe d line shows the 0.5 Jy 
limit of the AT20G Bright Source Sample ( Massardi et al. 20081) . 



be underestimated). For these sources, we find a mean flux ratio 

/ SaT20G \ 



\ SwMAP / 



1.01 it 0.03 . 



(7) 



The rms scatter of the individual fluxes is 0.31. The main contri- 
bution to the standard deviation in the flux ratios probably comes 
from variability, since the AT20G and WMAP measurements are 
generally not simultaneous. The typical uncertainty in the individ- 
ual WMAP and AT20G flux density measurements is 4 — 5 per cent. 
Since the formal standard error on the mean flux ratio is 0.028, we 
conclude that the AT20G and WMAP flux density scales are con- 
sistent to within 2-3 per cent at 20-25 GHz. 



8 ANALYSIS 
8.1 Polarisation 

To test the robustness of our polarization measurements we used 
the variability sample (see Section [231 , selecting a set of 142 ob- 
jects which had good quality observations in our 2004 October, 
2005 October and 2006 April observing epochs. Fig. Qj] shows the 
fractional polarisation at 20 GHz for the 2005 October versus 2004 
October observations. Crosses show measured polarisation values, 
and triangles show limits in one of the epochs. There is good corre- 
lation (albeit with significant scatter) between each pair of epochs. 

Fig. [T2l shows the 20 GHz polarised intensity versus the total 
flux density for variability sample sources observed in 2005 Octo- 
ber. The dashed horizontal line shows a 6 mJy limit in polarised 
flux, a level 3-4 times the typical error in polarised flux. Most of 
the limits (triangles) fall below this line, and so it was chosen as a 
reasonable level at which the polarised flux density can be reliably 
determined (see Section [33t . Although a small number of sources 
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Figure 13. Source density in equal area declination bands. Circles and 
crosses show sources above a 100 mJy and 50 mJy cutoff respectively. Er- 
rors are Poissonian. At a 50 mJy cutoff we have regions of incompleteness 
(notably 0°to -15° and -30° to -50°), as discussed in the text. 



Figure 11. Comparison of fractional polarisation at 20 GHz measured in our 
2004 October and 2005 October observing runs. Crosses show measured 
polarisation values, and triangles show limits in one of the epochs. 




\ 



20 GHz total flux density in 2005 (Jy) 

Figure 12. Polarised flux density versus total flux density at 20 GHz for the 
2005 October observing run. Crosses show measured polarisation values, 
and triangles show limits. 



have a measured polarised flux density below 6 mJy, further investi- 
gation found that these detections were near the 3a detection limit, 
and hence were not highly reliable. 

The triangle in the upper left of Fig. ITTI corresponds to source 
J084225-605350, which appeared highly polarized in 2005 (polar- 
ized flux 18 =b 1 mJy, fractional polarization 5.4 per cent) and 2006 
(16 =b 3 mJy, fractional polarisation 4.0 per cent) but not in 2004 
(< 6 mJy, fractional polarisation < 1.7 per cent). This deserves 



further investigation, as it raises the possibility that some AT20G 
sources may be genuinely variable in polarized flux. 

We found a good correlation between the fractional polarisa- 
tion measured at 5 GHz and 8 GHz, and also a reasonable corre- 
lation between the fractional polarisation measured at 5 GHz and 
20 GHz, although the large scatter in the latter case suggests there 
may be some change in the polarisation properties of the sources 
across the wider frequency range. There was also reasonable con- 
sistency in polarised position angles measured at each of the three 
frequencies. This will be analysed in more detail in Massardi et al. 
(in preparation). 



8.2 Source counts 

Fig.[l3]shows the 20 GHz source density with in equal area declina- 
tion bands, for flux density cutoffs of 100 mJy and 50 mJy. It shows 
that we have relatively homogeneous coverage of the southern sky 
down to 100 mJy. At a 50 mJy cutoff we have regions of incom- 
pleteness (notably 0°to —15° and —30° to —50°) as discussed in 
Section [6] This is primarily due to bad weather in the follow-up 
runs. 

Fig. Q3] shows a plot of the AT20G differential source counts 
(log(iV) — log (5)). These are well fit by a power law (f(x) = 
31x -2 ' 15 ) down to ~ 100 mJy where the flattening of the curve is 
a sign that the catalogue is becoming incomplete at this level. This 
agrees with our analysis in Section [6] Note that the source counts 
presented here are not corrected for completeness — this will be 
explored further in Massardi et al., (in preparation). 



8.3 Spectral index distribution 

In low-frequency radio surveys a commonly used diagnostic is the 
spectral index: 



log(Si/&; 
log(z/i/z/ 2 ) 



(8) 



where Si is the flux density measured at frequency vi and likewise 
for £2 and z/ 2 • 

The spectral index is only valid over frequency ranges in 
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Figure 14. log(iV) — log (S) source counts at 20 GHz, with Poisson errors. 
These are well fit by a power law down to ~ 100 mJy where the turnover 
indicates incompleteness. 



which the effects of spectral curvature can be ignored, which is 
clearly not the case for the AT20G sample. Hence for our analy- 
sis it is more appropr i ate to use a radio c olour-colour diagram (e.g. 
iKesteven et alJll977l ; [Sadler et al.ll2006h . Fig. [Bl shows a colour- 
colour plot for the 3763 AT20G sources with near- simultaneous 
data at 5, 8 and 20 GHz. The plot compares the lower frequency 
spectral index a (5, 8) with the higher frequency spectral index 
a(8,20). 

As noted previously bv lSadler et alj fcOOfj) and lMassardi et al.l 
(2008) this plot shows a wide range of spectral types at 20 GHz. 
Four main types can be identified: 

(i) Sources with steep (falling) spectra, which are shown in the 
lower left quadrant. These represent ~ 57 per cent of the sample 
and include a higher fraction of power law spectra. 

(ii) Sources with peaked (GPS) spectra, which rise at lower fre- 
quencies and fall at higher frequencies. These are shown in the 
lower right quadrant and represent ~ 21 per cent of the sample. 

(iii) Sources with inverted (rising) spectra, which are shown in 
the upper right quadrant. These make up ~ 14 per cent of sources. 

(iv) Sources with an upturn in their spectra, shown in the upper 
left quadrant. These represent ~ 8 per cent of the sample. 

The diagonal line indicates sources whose spectrum can be 
represented by a single power law. More sources lie below this line 
than above it, suggesting that most sources steepen with increasing 
frequency. For the flat spectrum sources in particular there is no evi- 
dence for power law spectra. This demonstrates that high frequency 
flux densities cannot be reliably estimated from low frequency sur- 
veys by extrapolating from a single power-law spectrum. 

We have also identified a new class of sources with ultra- 
inverted spectra (a(5, 20) > +0.7). These are shown as open cir- 
cles in Fig.[l5]and discussed in the next section. 



8.4 Low-frequency catalogues and Ultra-Inverted Spectrum 
(UIS) radio sources 

To test whether the AT20G survey has found any new radio-source 
population not seen at lower frequency, we cross-matche d the 
AT20G catalogue with the 1.4 GHz NV SS dCondon et alJl998h and 
843 MHz SUMSS jMauch et al J2003h and MGPS-2 jMurphv et alJ 
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Figure 15. Radio two-colour diagram. Crosses show AT20G sources ob- 
served at all three frequencies. Open circles show the Ultra-Inverted Spec- 
trum (UIS) sources with a spectral index of ok (5, 20) > +0.7. These 
sources are discussed in Section [84] 



l2007h surveys, which together cover the whole AT20G survey area. 
The 40 mJy cutoff of the AT20G catalogue is an order of magni- 
tude brighter than the NVSS and SUMSS/MGPS-2 completeness 
limits (which are typically 2.5-10 mJy, depending on declination), 
so we expect the vast majority of AT20G sources to have an NVSS 
or SUMSS match within 10-15 arcsec on the sky. AT20G sources 
without a counterpart in the NVSS and SUMSS catalogues must 
have either a sharply-rising radio spectrum (with a > +0.5) or 
highly variable radio emission. 

Our cross-matching identified 27 AT20G sources (0.4 per cent 
of the AT20G catalogue) with no counterpart in the NVSS and 
SUMSS/MGPS-2 catalogues, and these are listed in Table [7] For 
eleven of them, examination of the original NVSS and SUMSS sur- 
vey images shows a weak source at the AT20G position — these ob- 
jects are detected at the 3a level but fall below the NVSS/SUMSS 
catalogue limit. For the remaining sixteen objects, Table |7] lists a 
3a upper limit to the low-frequency NVSS/SUMSS flux density. 

Recent results from the 15 GHz 9C survey also confirm 
a population of high-frequency sources that fall below the lim- 
its of low-frequency catalogues. At flux densities of 10-15 mJy, 
Wald ramet ail 1 2009) find that 4.3 per cent of their 15 GHz sources 
are not listed in NVSS. 

All the sources in Table |7]have (non- simultaneous) 1-20 GHz 
spectral indices a(l, 20) > +0.7. An obvious question is whether 
these extreme spectral indices are the result of source variability 
between the NVSS/SUMSS and AT20G observing epochs. This 
does not generally appear to be the case. The sources in Table [7] 
are also among the most extreme AT20G objects in terms of their 
(simultaneous) 5-20 GHz spectral index, with a median value of 
a(5, 20) = +0.68 (the median a(5, 20) for the AT20G sample 
as a whole is -0.22). In all but two cases (J032945 -485420 and 
J123229-840247), the data in Tableware consistent with the ex- 
istence of a population of objects with a rapidly-rising power-law 
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Table 7. AT20G sources without a low-frequency counterpart in the SUMSS or NVSS catalogue. The (near-simultaneous) AT20G spectral index between 5 
and 20 GHz is also listed. Optical B(J) identifications and magnitudes are from the Supercosmos catalogue (Hambly et al. 2001), and infrared K magnitudes 
from the 2MASS Exten ded Source Catalog IJarrett et al.ll2000h . The redshift of J052755-471828 has been measured bv lDrake et all ( 120041) and the redshift 
of J09425 8 -604621 by Rad burn-Smith et al.N2006h . The median 20 GHz flux density of the sources in this table is 72 mJy. 



Flux density (mJy) 



AT20G Name 


S20 


S8 


S5 


NVSS 


SUMSS 


ex 


(1, 20) 


a(5,20) 


± 


Optical ID 


Comment 


J023611- 


-420337 


95 


75 


36 




7.3 




+0.81 


+0.68 


0.07 


Faint/blank 


No obvious optical ID 


J030406- 


-450342 


59 


49 


26 




< 3.3 


> 


+0.90 


+0.57 


0.09 


QSO? 


Stellar ID, B(J)=18.0 mag 


J032945- 


-485420 


40 


58 


47 




< 3.3 


> 


+0.79 


-0.11 


0.08 


QSO? 


Stellar ID, B(J)=19.1 mag 


J034258- 


-431813 


97 


97 


39 




9.4 




+0.74 


+0.64 


0.09 


QSO? 


Stellar ID, B(J)=1 8.4 mag 


J044023- 


-473218 


61 


21 


9 




< 3.3 


> 


+0.92 


+1.34 


0.12 


Faint/blank 


Possible faint ID, B(J)=22.7 mag 


J052755- 


-471828 


106 








9.9 




+0.75 






Galaxy 


2MASS galaxy K=13.0mag, z=0.134 


J054417- 


-641914 


123 


66 


32 




< 3.3 


> 


+1.14 


+0.94 


0.08 


Faint/blank 


Crowded field near LMC, no obvious ID 


J070949- 


-381152 


86 


33 


6 


< 2.5 


< 5.1 


> 


+1.33 


+1.87 


0.15 


Galaxy 


2MASS galaxy with K=12.7mag 


J073040- 


-544152 


70 








< 7.5 


> 


+0.70 






Faint/blank 


Low S/N in SUMSS image, no obvious ID 


J080931- 


-472011 


86 


50 


27 




7.7 




+0.76 


+0.81 


0.09 




Low Galactic latitude, b=-7.7° 


J094258- 


-604621 


51 


33 


23 




8.7 




+0.56 


+0.59 


0.06 


Galaxy 


2MASS galaxy K=12.9mag, z=0.096 


J095159- 


-183703 


70 


71 


38 


< 3.3 




> 


+1.15 


+0.43 


0.08 


Faint/blank 


No obvious optical ID 


J104227- 


-210556 


58 






< 2.0 




> 


+1.27 






Faint/blank 


No obvious optical ID 


J111015- 


-665531 


136 


66 


21 




6.0 




+0.99 


+1.31 


0.07 


QSO? 


Stellar ID B(J)=19.1 mag 


J111246- 


-203932 


78 


30 


6 


< 2.0 




> 


+1.38 


+1.80 


0.17 


Faint/blank 


No obvious optical ID 


Jl 11605- 


-263758 


40 


5 


3 


< 1.2 




> 


+1.32 


+1.82 


0.34 


Faint/blank 


No obvious optical ID 


Jl 14844- 


-781933 


124 


63 


22 




< 5.0 


> 


+1.01 


+1.21 


0.07 


Galaxy 


B(J)=18.5 mag galaxy 


J123229- 


-840247 


48 


86 


79 




6.3 




+0.64 


-0.35 


0.08 


Faint/blank 


No obvious optical ID 


J140257- 


-664031 


72 


58 


26 




< 5.0 


> 


+0.84 


+0.71 


0.07 




Low Galactic latitude, b=-4.8° 


J161845- 


-142428 


84 






< 2.5 




> 


+1.32 






Faint/blank 


No obvious optical ID 


J171043- 


-471820 


76 








6.8 




+0.76 








Low Galactic latitude, b=-4.4° 


J173039- 


-211112 


54 






< 2.5 




> 


+1.16 






Faint/blank 


No obvious optical ID 


J181225- 


-712006 


43 


39 


26 




5.2 




+0.67 


+0.35 


0.09 


Galaxy 


2MASS galaxy K=12.9mag 


J195906- 


-755202 


42 


30 


23 




< 3.5 


> 


+0.78 


+0.42 


0.07 


QSO? 


Stellar ID, B(J)=20.6mag 


J200012- 


-474951 


78 


64 


28 




7.4 




+0.74 


+0.72 


0.10 


QSO? 


Stellar ID, B(J)=20.2 mag 


J205503- 


-635207 


44 


29 


12 




4.6 




+0.71 


+0.91 


0.09 


QSO? 


Stellar ID, B(J)=20.2 mag 


J220413- 


-465424 


99 


103 


30 




< 4.0 


> 


+1.01 


+0.84 


0.08 


Faint/blank 


No obvious optical ID 



radio spectrum between 1 and 20 GHz and little or no variability 
on timescales of up to a decade. 

By analogy with the Ultra-Steep spectrum (USS) radio 
sources with a < — 1.3 felumenthal & Milevlll979l :l Tielens et al.1 
Il979h we introduce the term "Ultra-Inverted Spectrum (UIS) radio 
source" to describe the class of radio sources with a spectral index 
a(5,20) > +0.7. 

The AT20G catalogue contains 45 sources with a (5, 20) > 
+0.7 (roughly 1.2% of the AT20G sources with 5 and 8 GHz data), 
and these are listed in Table [8] A few of these sources are detected 
in linear polarization, confirming that the emission mechanism is 
non-thermal, so we assume for now that all the AT20G UIS ob- 
jects are extragalactic non-thermal sources. Further investigation is 
needed to confirm this, but we note that the UIS sources have the 
same distribution in Galactic latitude as the main AT20G sample 
and show no concentration towards the Galactic plane. 

Fig. [T6l compares the AT20G sources in Table [7]with the sam- 
ple of High Frequency Peakers (HFPs) studied by Dallacasa et al. 
(2000), who measured simultaneous flux densities at several fre- 
quencies between 1.4 and 22.5 GHz for a sample of 55 sources 
selected to have inverted spectra with a(l, 5) > +0.5. Our AT20G 
UIS sources with a(5, 20) > +0.7 have radio spectra which 
rise more rapidly with frequency than any of the HFP sample. If 
the AT20G UI S sour c es pea k at frequencies near or above above 
20 GHz, the lO'Deal Jl998h relation between spectral peak and 
source size implies that these are expected to be very young, com- 
pact radio sources less than a few tens of parsecs in size. 

Only a handful of similarly-extreme inverted- spectrum radio 
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Figure 16. Comparison of 1 — 5 and 5 — 20 GHz spectral indices for the 
so urces in T able [7] (filled circles) and the 'High-Frequency Peaker' sample 
of lDallacasa et al J 12000 ) (crosses). Open circles show AT20G UIS sources 
(Table[8j which have both a(5, 20) > +0.7 and a(l, 20) > +0.7. 
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Table 8. "Ultra-Inverted Spect rum" (UIS) AT20G sources with 5-20 GHz spectral index a ^ +0.7. 'He07' in the Notes column indicates a source which is 
also in the CRATES sample of Healey et al. (2007). Identifications an d redshifts are fro m the NASA Extragalactic Database (NED). Redshift references are: 
Su04 = ISulentic et"ai1 120041) ; St91 = IStocke et alj fl99lh : 6dFGS = I Jones etal] l2009h : 2dFGRS = |Colless etal.ll2 001). The strongest source in this table, 
J132649-525623 (=PMN J1326-5256) has been observed to show intra-day variability at 6.6 GHz iMcCulloch et al.ll2005l) . 
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sources have previousl y been identified. They include the z = 
0.089 galaxy III Zw 2 (iFalcke et alJll999fc iBrunthaler et al.ll2003l) 
and t h e z 0.644 quasar RXJ 1415+3337 (lOrienti & Dallacasal 
2008). IFalcke et all Jl999h measured a spectral index of a = 
+ 1.9 ±0.1 between 5 and 10 GHz for III Zw 2, with a spectral peak 
near 40 GHz, while RXJ 1415+3337 has a = +0.9 between 1 and 
22 GHz. Both sources show an evolution of the radio spectrum over 
timescales of several years, with the spectral peak slowly moving 
to lower frequencies, and III Zw 2 can be modelled in terms of 
a single, adiabatically-ex panding homogeneous radio component 
(lOrienti & Dallacasall2008b . 



8.5 Optical Identifications 

Optical counterparts of the full AT20G sample were found by 
cross-matching the radio positions with o ptical positions in the Su- 
perCOSMOS database dHamblv et al.l200ll) . AT20G objects within 
10° of the Galactic plane were excluded from this analysis due to 
the presence of foreground stars and Galactic dust extinction, leav- 
ing a total of 4932 objects used in this analysis. Optical identifi- 
cations were chosen to be the closest optical source to the radio 
position within a 2.5 arcsec radius and brighter than a B magnitude 
of 22, which is the SuperCOSMOS completeness limit. Previous 
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studies have found that 97 per cent o f sources selected in this way 
are likely to be genuine associations ( Sadler et al. 20q3). 

Using this selection method there are 2958 AT20G sources 
(60 per cent) with optical identifications, the majority of which 
are QSOs. This is much higher than seen in low frequency ra- 
dio surveys which typically have an optical identification rate of 
25—30 per cent. This again highlights the significant difference of 
the AT20G source population from other radio selected AGN sam- 
ples. Full analysis of the optical properties of these objects will be 
published separately (Mahony et al. in preparation) 



9 CONCLUSIONS AND FUTURE WORK 

We present a catalogue of 5890 sources from the Australia Tele- 
scope 20 GHz survey, the deepest large-area survey at high radio 
frequency. For 3766 of these sources we have near- simultaneous 
5 and 8 GHz measurements, and 1559 sources have a detection in 
polarized total intensity at one or more of the three frequencies. 

The 20 GHz flux densities measured for the strongest AT20G 
sources are in excellent agreem ent with the WMAP 5 -year source 
catalogue recently published by I Wright et al.l (120091) . and we find 
that the WMAP source catalogue is close to complete (and highly 
reliable) for sources stronger than 1.5 Jy at 23 GHz. 

We identify a population of Ultra-Inverted Spectrum radio 
sources with a spectral index of a (5, 20) > +0.7. These are rare 
sources, comprising roughly 1.2 per cent of the AT20G population. 

There are several ongoing projects as part of the AT20G. Mas- 
sardi et al. (in preparation) will present a statistical analysis of the 
AT20G sources, and Hancock et al. (in preparation) will present 
results from the original scanning survey. Mahony et al. (in prepa- 
ration) is carryi ng out an analysis o f the optical identifications of 
AT20G sources. I Sadler et all (12008b are following up s ubsamples 
of the AT20G sources at 95 GHz. iMurphv et all 12009) is target- 
ing a subset of Galactic sources that were excluded from the main 
follow-up survey. Chhetri et al. (in preparation) has carried out a 
search for gravitational lens candidates and planetary nebulae using 
data from the 6 km baseline. In a related project, Sadler et al. (in 
preparation) is conducting a deeper survey at 20 GHz to explore the 
high frequency radio population at much lower flux densities. 
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APPENDIX A: WMAP SOURCES MISSING FROM THE 
AT20G CATALOGUE 

Table IA1I lists the eight WMAP sources that do not have a cata- 
logued bright (> 250 mJy) AT20G source within 21.35 arcmin of 
the WMAP position. These 'missing' sources fall into several cat- 
egories which are discussed in the following sections. 



Al Bright nearby radio galaxies with large angular size 

As noted by Massar di et al ] d2008h . the radio galaxy Fornax A 
(WMAP J0322-3711) is not detected in the AT20G survey be- 
cause most of its 20 GHz flux density arises from diffuse emission 
associated with the lobes, which are resolved out by the 2 arcmin 
ATCA beam. Fornax A is the only bright extragalactic source 
known to be missing from the AT20G catalogue. 

Two other WMAP sources (J0133-3627 and J0636-2031) 
correspond to peaks in the extended emission of the nearby radio 
galaxies NGC612 and PKS 0634-20 respectively. Both these 
galaxies are detected by AT20G, and the large WMAP-AT20G 
position offsets arise from the complex structure of the extended 
radio emission. 



A2 Extended Galactic sources 

Two WMAP objects (J0519-0539 and J2 229- 2050) a re as - 
sociated with extended Galactic sources. Wrig ht et all |2009) 
note that J0519-0539 is a blend of two Lynds Bright Nebulae, 
LBN207.65-23.il and LBN 207.29-22.66, while J0636-2031 
is part of the Helix Nebula, NGC7293. Both objects are signif- 
icantly larger than the AT20G beam, and appear to be mostly 
resolved out in our survey. 



A3 Extragalactic WMAP sources not found in AT20G 

There are three remaining WMAP sources (J0632— 6928, 
Jl 150— 7927 and J1637— 7714) for which no obvious counterpart 
can be found in the AT20G catalogue. All three of these sources 
are detected in each of the five WMAP single-year images (Wright 
et al. 2009), and do not appear to vary significantly in flux density 
over this five-year timespan. 

In two cases (WMAP J0632-6928 and Jl 150-7927), there 
is no catalogued AT20G source within the 0.9° WMAP beam at 
the listed WMAP position. Unless these WMAP sources are very 
diffuse, it seems unlikely that they are real. 

J0632-6928 has a WMAP flux density of 0.4 Jy, making it 
the weakest K-band (23 GHz) detection in the Wright et al. (2009) 
catalogue, but it is detected (as a sub-Jy source) in all 5 WMAP 
bands. The closest AT20G source, J063455-694532, is 0.4 deg 
away and has a flux density of 66mJybeam _1 at 20 GHz and 
154 mJy beam -1 at 5 GHz. Its PMN counterpart is PMN J0634- 
6945 with a flux density of 92±8 mJy beam" 1 , and the SUMSS 
flux density is 123. 1 ±3.8 mJy beam -1 . AT20G J063455-694532 
seems too faint to be a plausible identification for the WMAP 
source. 

The third sour c e (W MAP J1637-7714) is identi- 
fied by Wrig ht et all (2009) with a nearby 5 GHz source 
PMN J1636— 7713, which is also detected as a weak 
(40±2mJybeam" 1 at 20 GHz) AT20G source. Once again, 
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however, the catalogued AT20G flux density is well be- 
low the WMAP detection limit. A stronger AT20G source 
(AT20G J164416-771548, with a 20 GHz flux density of 
399±20mJybeam _1 ) lies 22arcmin from the WMAP position, 
and may be responsible for the WMAP detection. 
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Table Al. The eight southern WMAP sources that are not matched with a bright (> 25 mJy beam 1 ) AT 20G source, within a matching radius of 
21.35 arcmin. The listed WMAP K-band (23 GHz) flux densities are from the 5-year catalogue (Wright etal. 2(X)3). 



WMAP position 


WMAP 


± 


AT20G 


Notes 


(J2000) 


s 23 (Jy) 


(Jy) 


s 20 (Jy) 




01 33 26.6 -36 27 09 


0.6 


0.1 




WMAP source is a hotspot of the radio galaxy NGC 612, see text 


03 22 25.4 -37 11 25 


18.5 


3.1 




Radio galaxy Fornax A, see text 


05 19 21.6 -05 39 37 


2.4 


0.1 




Blend of Galactic emission nebulae (Wright et al. 2009), no AT20G counterpart 


06 32 21.1 -69 28 32 


0.4 


0.0 




No AT20G counterpart 


06 36 31.8 -20 31 38 


1.1 


0.0 


0.13 


Part of the radio galaxy PKS 0634—20, matched with two faint AT20G sources 


11 50 12.9 -79 27 35 


1.2 


0.0 




No AT20G counterpart 


16 37 52.5 -77 14 58 


1.4 


0.1 


0.04 


Faint AT20G counterpart, see text 


22 29 47.1 -20 50 28 


0.9 


0.1 




Part of the Galactic planetary nebula NGC 7293, no AT20G counterpart 



Table 7. The first 30 sources in the AT20G source catalogue. 



IAU designation 


a 


8 


S20 


^20 


s 8 




s 5 


^5 


Ep a 


Q b 


F c 


P20 <y p 20 


m 


PA 20 


Ps crp 8 


m 


PA 8 


P 5 




m 


PA 5 






(J2000) 


(J2000) 


(mJy) 


(mJy) 


(mJy) 








(mix) 




(%) 


(°) 


(mJy) 


(%) 


(°) 


(mJy) 


(%) 


(°) 


AT20GJ000012- 


-853919 


00:00:12.78 


-85:39:19.9 


98 


5 


63 


4 


63 


4 


222 


g 




<10 


- 


<9.8 




<6 - 


<9.6 




<6 




<9.5 




AT20GJ000020- 


-322101 


A A aa r\r\ o o 

00:00:20.38 


O O O 1 A 1 O 

—32:21:01.2 


118 


6 


315 


16 


515 


26 


111 


g 




<6 


- 


<5.1 




<6 - 


<1.9 




<6 




<1.2 




AT20GJ000105- 


-155107 


00:01:05.42 


-15:51:07.2 


297 


19 


295 


15 


257 


13 


444 


g 




<8 


- 


<2.7 




<6 - 


<2.0 




<6 




<2.3 




AT20GJ000106- 


-174126 


00:01:06.31 


— 17:41:26.2 


73 


8 










6 • • 


g 




<6 


- 


<8.6 


















AT20GJ000118- 


-074626 


A A r\ i i o f\ A 

00:01:18.04 


AO /I O O 

—07:46:26.8 


177 


9 










6 • • 


g 




<8 


- 


<4.4 


















AT20GJ000124- 


-043759 


00:01:24.50 


-04:37:59.6 


50 


3 










6 • • 


g 




<6 


- 


<12.6 


















AT20GJ000125- 


-065624 


aa ai r\ c c r\ 

00:01:25.59 


—06:56:24.7 


77 


4 










6 • • 


g 




<8 


- 


<10.2 


















AT20GJ000212- 


-215309 


A A AO i o AO 

00:02:12.02 


O 1 C O A A A 

—21:53:09.9 


165 


11 










4 • • 


g 




<6 


- 


<3.6 


















AT20GJ000221- 


-140643 


A A AO o 1 ni 

00:02:21.71 


1 /I A A O /A 

— 14:06:43.9 


48 


3 










7 • • 


g 




<6 


- 


<12.4 


















AT20GJ000230- 


-033140 


A A AO o a s a 

00:02:30.60 


AO O 1 /I A 1 

—03:31:40.1 


53 


3 










6 • • 


g 




<9 


- 


<17.5 


















AT20GJ000249- 


-211419 


00:02:49.85 


O 1 1 /l 1 A O 

—21:14:19.2 


100 


7 










4 • • 


g 




<7 


- 


<6.9 


















AT20GJ000252- 


-594814 


A A AO c o a o 

00:02:52.93 


r- A /in 1 A f\ 

—59:48:14.0 


71 


3 


64 


3 


57 


3 


222 


g 




<9 


- 


<12.2 




<6 - 


<9.4 




<6 




<10.5 




AT20GJ000253- 


-562110 


A A AO c o C 

00:02:53.65 


C f O 1 1 A O 

—56:21:10.8 


94 


5 


229 


12 


403 


20 


222 


g 




<8 


- 


<8.3 




<6 - 


<2.6 




<6 




<1.5 




AT20GJ000303- 


-553007 


A A AO r\ O A C 

00:03:03.45 


r" r" o A AO 1 

—55:30:07.1 


44 


3 


48 


3 


52 


3 


222 


g 




<7 


- 


<17.0 




<6 - 


<12.6 




<6 




<11.5 




AT20GJ000311- 


-544516 


r\r\ AO -i -i r\ a 

00:03:11.04 


—54:45:16.8 


95 


3 


313 


3 


552 


9 


222 


g 


e 


<8 


- 


<19.9 




<6 - 


<11.8 




15 


2 


8.0 


—28 


AT20GJ000313- 


-590547 


r\r\ AO 1 o o o 

00:03:13.33 


r" A A r" /I o o 

—59:05:47.7 


49 


3 


101 


5 


151 


8 


222 


g 




<8 


- 


<17.2 




<6 - 


<5.9 




<6 




<4.0 




AT20GJ000316- 


-194150 


A a AO 1 A 

00:03:16.06 


1 A A 1 r" A O 

— 19:41:50.7 


76 


5 


162 


9 


182 


9 


444 


g 




<8 


- 


<10.6 




<7 - 


<4.4 




<6 




<3.3 




AT20GJ000322- 


-172711 


A A AO OO A r" 

00:03:22.05 


1 O O O -I -I A 

— 17:27:11.9 


386 


43 










6 • • 


g 




14 


2 


3.7 


18 
















AT20GJ000327- 


-154705 


A A AO OH O C 

00:03:27.35 


1 r" /in rvr A 

— 15:47:05.4 


129 


8 










4 • • 


g 




<10 


- 


<7.5 


















AT20GJ000404- 


-114858 


A A A A A A O O 

00:04:04.88 


1 1 /in r- o A 

— 11:48:58.0 


680 


32 










6 • • 


g 




<7 


- 


<1.0 


















AT20GJ000407- 


-434510 


A A A /I AO O A 

00:04:07.24 


/I O /I C 1 A A 

—43:45:10.0 


199 


10 


211 


11 


244 


12 


111 


g 




<8 




<4.1 




<6 - 


<2.8 




<6 




<2.5 




AT20GJ000413- 


-525458 


00:04:13.97 


-52:54:58.7 


65 


3 


98 


4 


192 


6 


222 


g 


e 


<7 




<15.7 




<6 - 


<14.2 




<6 




<8.4 




AT20GJ000435- 


-473619 


00:04:35.65 


-47:36:19.0 


868 


36 


970 


49 


900 


45 


111 


g 




15 


3 


1.7 


-51 


30 2 


3.1 


-45 


25 


1 


2.8 


-43 


AT20GJ000505- 


-344549 


00:05:05.94 


-34:45:49.6 


131 


6 


142 


7 


134 


7 


111 


g 




<6 




<4.6 




<6 - 


<4.2 




<6 




<4.5 




AT20GJ000507- 


-013244 


00:05:07.03 


-01:32:44.6 


81 


4 










6- 


g 




<6 




<7.8 


















AT20GJ000518- 


-164804 


00:05:18.01 


-16:48:04.9 


142 


9 










4.. 


g 




<6 




<4.2 


















AT20GJ000558- 


-562828 


00:05:58.32 


-56:28:28.9 


151 


5 


376 


5 


677 


11 


222 


g 


e 


<7 




<7.2 




<6 - 


<8.5 




<6 




<2.7 




AT20GJ000600- 


-313215 


00:06:00.47 


-31:32:15.0 


63 


4 


53 


3 


52 


3 


111 


g 




<9 




<13.8 




<6 - 


<11.3 




<6 




<11.6 




AT20GJ000601- 


-295549 


00:06:01.14 


-29:55:49.6 


97 


6 


187 


10 


228 


10 


444 


g 




<6 




<6.2 




<7 - 


<3.9 




<6 




<2.6 




AT20GJ000601 


-423439 


00:06:01.95 


-42:34:39.8 


110 


5 


259 


14 


532 


27 


111 


g 




11 


3 


9.8 


32 


14 1 


5.5 


22 


<6 




<1.1 





NOTES: 

a Epoch of follow-up observations for 20 GHz, 8 GHz and 5 GHz fluxes respectively. The epochs are listed in Table 3. 
b Quality flag: g (good) or p (poor). c Other flags: see Section[5]for description. 



